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Abstract. 

The Sun's gradual brightening will seriously compromise the Earth's biosphere 
within ~ 10^ years. If Earth's orbit migrates outward, however, the biosphere could 
remain intact over the entire main-sequence lifetime of the Sun. In this paper, we 
explore the feasibility of engineering such a migration over a long time period. The 
basic mechanism uses gravitational assists to (in effect) transfer orbital energy from 
Jupiter to the Earth, and thereby enlarges the orbital radius of Earth. This transfer 
is accomplished by a suitable intermediate body, either a Kuiper Belt object or a 
main belt asteroid. The object first encounters Earth during an inward pass on its 
initial highly elliptical orbit of large (~ 300 AU) semimajor axis. The encounter 
transfers energy from the object to the Earth in standard gravity- assist fashion by 
passing close to the leading limb of the planet. The resulting outbound trajectory 
of the object must cross the orbit of Jupiter; with proper timing, the outbound 
object encounters Jupiter and picks up the energy it lost to Earth. With small 
corrections to the trajectory, or additional planetary encounters (e.g., with Saturn), 
the object can repeat this process over many encounters. To maintain its present 
flux of solar energy, the Earth must experience roughly one encounter every 6000 
years (for an object mass of 10^^ g). We develop the details of this scheme and 
discuss its ramifications. 

Keywords: Orbits, Celestial mechanics 



1. Introduction 



As the Sun burns through its hydrogen on the main sequence, it steadily 
grows hotter, larger, and more luminous. Stellar evolution calculations 
show that in -^1.1 billion years the Sun will be 11% brighter than it 
is today (e.g., Sackmann et al., 1993). Global climate models indicate 
that such an increase in insolation would drive a "moist greenhouse" 
on the Earth (Kasting, 1988; Nakajima et al., 1992) which will have 
a catastrophic effect on the surface biosphere. In 3.5 billion years, the 
total luminosity of the Sun will be 40% larger than the present value. 
Under such conditions, the Earth will undergo a catastrophic "runaway 
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greenhouse" effect (Kasting, 1988), which will likely spell a definite end 
to life on our planet. 

Although the Earth's ecosystem will be seriously compromised within 
a billion years, the Sun is presently less than halfway through its main 
sequence life. Indeed, in 6.3 billion years, the luminosity of the Sun is 
expected to be "only" a factor of 2.2 greater than its current value. At 
that time, a planet located at 1.5 AU from the Sun would receive the 
same flux of solar energy that is now intercepted by the Earth. 

If the radius of the Earth's orbit were somehow to be gradually in- 
creased, catastrophic global warming could be avoided, and the lifespan 
of the surface biosphere could be extended by up to five billion years. 
In this paper, we study the feasibility of altering planetary orbits over 
long time scales. Special attention will be given to the specific case 
of the Earth, but many of the issues we address are of more general 
astronomical and astrobiological interest. 

The present orbital energy of the Earth is —2.7 x 10^^ erg. Moving 
the Earth to a circular orbit of 1.5 AU radius would require 8.7 x 
■|^q39 gj,g ^j-^ attractive scenario for gradually increasing the Earth's 
orbital radius is to successively deflect a large object or objects from 
the outer regions of the solar system (the Oort Cloud or the Kiiiper 
Belt) onto trajectories which pass close to the Earth. By analogy to the 
gravity-assisted flight paths employed by spacecraft directed to outer 
solar-system targets (e.g.. Bond and Anson, 1972, Minovitch, 1994), the 
close passage of such an object to the Earth can result in a decrease 
in the orbital energy of the object and a concomitant increase of the 
Earth's orbital energy. For optimal trajectories which nearly graze the 
Earth's atmosphere, the energy boost imparted to the Earth is 2.4 x 10^^ 
erg gm~^ of object mass (Niehoff, 1966). Work by Sridhar and Tremaine 
(1992) suggests that even bodies that are weakly held together ("rubble 
piles" ) can survive passages that approach less than 1 Earth radius from 
the Earth's surface, allowing energy transfers of 10^^ erg gm~^. 

Typical masses for large Kuiper belt objects are of the order of 10^^ 
grams, meaning that roughly 10^ passages (involving a cumulative flyby 
mass of approximately 1.5 Earth masses) would be needed to move the 
Earth out to 1.5 AU. Thus, over the remaining lifespan of the Sun, 
approximately one passage every 6000 years on the average would be 
required. 

The outer reaches of the Solar System contain an ideal reservoir of 
material which could be used to move the Earth. The Kuiper Belt is 
populated by a large number of objects that are larger than 100 km 
in diameter; the Kuiper belt may contain as many as 10^ such bodies, 
totaling perhaps 10% of the Earth's mass (Jewitt, 1999), although these 
numbers remain uncertain. The Oort cloud is believed to contain about 
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10^^ objects totaling 30 or more Earth masses (see, e.g., Weissman, 
1994). As evidenced by the frequent passage of long period Sun-grazing 
comets originating in this region, many Oort cloud objects would need 
only small trajectory changes in order to bring them into appropriate 
Earth-crossing orbits. Indeed, strategies for modifying the orbits of 
asteroids and comets have been extensively discussed in the context 
of mitigating the hazard posed by such objects impacting the Earth 
(see, e.g., Ahrens and Harris, 1992, Melosh et al., 1994, Solem, 1991). 
Alternatively, a main belt object could be deflected into an orbit which 
has an aphelion in the outer solar system. 

Our approach in this paper is as follows: In §2, we discuss the details 
of our gravity assist scheme. This scheme uses an asteroid or large 
comet as a catalyst to transfer orbital angular momentum and energy 
from Jupiter to the Earth. Wc investigate the energy requirements 
of the scheme, the nature of the course corrections demanded, and 
also the needed accuracy. In §3, we discuss additional considerations, 
such as long term orbital stability, complications produced by other 
planets, and larger issues. We present our conclusions in §4. Although 
this problem raises many possible interesting (and rather speculative) 
issues, the present paper discusses only a few of them. 



2. The gravity-assist scheme 

As mentioned in the introduction, our underlying scenario uses re- 
peated gravity assists to (in effect) transfer orbital energy from Jupiter 
to the Earth, thus enlarging the Earth's orbit and reducing the received 
solar flux. Multiplanet encounter trajectories have been discussed for 
more than 25 years (e.g.. Bond and Anson, 1972) and are now com- 
monplace features of interplanetary exploration, as evidenced by the 
Galileo and Cassini missions. 

The underlying dynamics of the scheme are shown in Figs. 1 and 2. 
The object "O", a suitable Kuiper belt object or main belt asteroid, 
flrst encounters the Earth during an inward pass on its initial highly 
elliptical orbit of large (9(300 AU) semimajor axis. The encounter trans- 
fers energy from O to the Earth in standard gravity-assist fashion by 
passing close by the leading limb of the Earth. The resulting orbit of O 
then crosses the orbit of Jupiter; with proper timing, it will encounter 
Jupiter on its outbound swing, pass by Jupiter, and regain the energy it 
lost to the Earth. It would appear, however, that O also gains angular 
momentum relative to its incoming orbit, so that the return orbit is 
less elliptical than the initial one. As a result, a modest amount of 
energy must be expended to restore O's orbit to its initial parameters. 
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Figure 1. Geometry of orbital encounters and gravity assists for a generic encounter 
of the O and the Earth, a) Encounter in the heliocentric frame. The distance scale 
is in AU. b) Encounter in the planet-centered frame, in which the object's flyby past 
the planet is a hyperbola with close-approach distance h and turning angle a. The 
distance scale is in units of 10^ cm, for a typical Earth encounter. 



unless further encounters are included that can minimize the required 
expenditure. 

As discussed below, larger orbits for O entail lower first-order energy 

requirements. On the other hand, an orbit that is too large will have 
too long a period to work well in the scheme (unless multiple objects 
are used). For purposes of discussion, we use an orbit with a semi- 
major axis of 325 AU, whose period of 5859 years is compatible with 
the ~ 6000 year average period between encounters. In fact, the Sun's 
brightening is slow at first and then speeds up; thus the optimum orbit 
size will decrease with time. 

2.1. Formulation 

We begin by assuming that Earth and Jupiter are on circular orbits of 
zero inclination with radii of 1 and 5.2 AU, respectively. Obviously, 
a more detailed study would have to take into account the actual 
elements of the planetary orbits. To outline the scheme, however, the 
idealized case is adequate. The calculation of energy transfer and orbital 
parameters is made using the so-called "patched-conic" approximation 
(Battin, 1987, Bond and Allman, 1986). In this approximation, the 
orbit of the object is treated as a series of two-body problems. Far 
from planets, O follows a Keplerian ellipse about the Sun. At plane- 
tary encounters, O's path is given by a two-body scattering encounter. 
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Figure 2. Layout of successive encounters for the Earth- Jupiter scheme, for an 
orbit with initial apheliou Ro at 650 AU and aphelion tangential velocity Vo ~ 6000 
cm s~^. Note the changes of scale from frame to frame, a) Initial orbit to Earth 
encounter b) Orbit post-Earth to Jupiter c) Orbit post-Jupiter d) Initial (dotted) 
and return (solid) orbits compared. 



Numerical integrations of the full four body system indicate that this 
approximation gives perfectly adequate results. 

The incoming orbit of O (characterized by subscripts 0) is param- 
eterized for convenience by an aphelion distance Rq and tangential 
velocity Vq. Then the incoming angular momentum and energy are 
given by 



Co = RoVq, 



(1) 
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where fiQ = GMq. Assuming ho < 0, the semi-major axis, eccentricity, 
and longitude of perihehon of the orbit are 

1/2 



— , eo - i H 2~ 

2^0 V f^e 



«o = -^, eo=(l + — ^1 , and a;o = 0. (2) 



For an Earth encounter, the eccentricity must be large enough so that 

the perihelion distance ao(l — eo) < i?©, the orbital radius of the 
Earth. If this constraint is satisfied, the Earth encounter takes place at 
longitude (pE, where 

cos = — - and po = — ■ (3) 



eo V-Re / Mo 
The detailed geometry of an encounter is illustrated in Fig. 1. For 

an "incoming" (pre-pcrihclion) encounter, it < (p^ < 2tt. The Earth's 
orbital velocity is V® = {^q/ R^)^/"^ . The object's speed Vg, and tan- 
gential and radial velocities Vte,Vre at the encounter follow from 
conservation of angular momentum and energy: 

1/2 

/r> _ T^2 T 

^TE- 



Ve 



Vte = co/i?®, V^e = Vi - 



(4) 

In the Earth's frame of reference, the velocities are V^^ = Vte — 

V®, Vl^E = Vre, and the encounter speed is given by (V^)^ = (V^y^;)^ + 
(y^^)^. The velocity vector of the encounter in the Earth's frame makes 
an angle Pe with respect to the orbital velocity of the Earth, where 

cosPe = V^e/Ve- 

In the two-body treatment, the effect of the encounter is to turn the 

velocity vector (in the Earth frame) through an angle as, where aE 

depends on the encounter velocity and the impact parameter Be- The 

encounter can be timed to produce a specified minimum distance from 

the Earth, bE- The minimum distance, the impact parameter Be, and 

the turning angle aE are related by 

(5) 

For an encounter in which O loses energy and Earth gains energy, 
aE > 0. The post-encounter tangential velocity, from which the energy 
transfer is found, is then given by V^^e = V^cos(/3e + ag). Similarly, 
the post-encounter radial velocity is V^^e = ^e sin(/3£; + aE)- 

The change in energy per unit mass of the object, from pre-encounter 
to post-encounter, is then given by 

AQb = (1/2)[(V£ • V£)post - (Vb • V£)pre] = V^{V^E ' V^e)- (6) 
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The corresponding change in the Earth's orbital energy is thus —MqAQe, 
where Mq is the mass of O. As mentioned above, AQe will be neg- 
ative (and hence the Earth will gain energy) if O passes "in front" 
of the Earth. The amount of energy transfer depends not only on the 
minimum approach distance but also on the encounter geometry, i.e., 
Pe, and the encounter speed V^, which in turn depend on the longitude 
(j)E of the encounter. Generally speaking, the most effective encounters 
occur for near but not quite "grazing" encounters for which ~ ±0.5 
rad, not far from O's perihelion. If bE can be taken as small as 10^ cm 
(about 1.6 Earth radii), encounter transfer energies of up AQe 10^^ 
erg gm~^ can be achieved. This value is approximately 60% of the 
maximum AQe = V"0V^irc, where V^irc is the circular velocity in Earth 
orbit at a radius Be from the Earth's center. 

We denote post-Earth-encounter quantities by the subscript 1. For 
these post-Earth variables, we use Cartesian vectors in the orbital 
plane, for which = {Xq,Y^) = {R^ cos (I)e,R® sin (f)E)- Similar 
formulae obtain for the post-Earth velocity (in the solar frame) Vi. In 
particular, 

Vxi = V^E cos (l)E-{VTE+y®) sin (^E, Vyi = V^E sin ^e+{Vte+V(s) cos (I)e. 

(7) 

The angular momentum is then given by Ciz = Rgj x Vi, and the 

Laplace vector Pi = — jl(.^YiQ / R(q + ciz x Vi. With these forms, we 
obtain the orbital elements of the object O in its post-Earth orbit, i.e., 

hi = — — — , ai = -TTT-, ei = , wi = tan {Pyi/Pxi) 

I zrii jiQ 

(8) 

Examination of the post-Earth orbital elements shows that the new 
orbit of O crosses the orbit of Jupiter. We may therefore schedule an 
encounter with Jupiter to regain energy lost by O to Earth. As before, 
treating the orbit of Jupiter as a circle and the orbit of O as an ellipse 
in the plane, we find that the longitude of the Jupiter encounter is 
given by 



cos((/)j -u)i) = — — 1 , where pi = — . (9) 

ei \R% ) He 

The encounter of O with Jupiter implies similar considerations (with 

respect to the change in orbital parameters) as the encounter with 
Earth. The orbital elements give us the tangential encounter velocity 
Vj^j (in Jupiter's frame) and the encounter speed Vj. Our initial idea 
was to set the encounter geometry so as to yield an energy gain AQj 
(by O) equaling the amount lost at the Earth. As noted below, however. 
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a more efficient encounter (in terms of tfie final velocity change of O) 
is one that yields a post- Jupiter orbit with an aphelion equal to the 
original value Rq. In that case, it is simplest to search numerically 
for the desired Jupiter encounter distance bj. The encounter geometry 
gives us Pj; the impact parameter and aj then follow from bj as above. 

Finally, we must work out the orbital elements OR,eR,u;R for O on 
its post-Jupiter return orbit. The succession of encounters is shown in 
Fig. 2, for our example orbit with aphelion at 650 AU and and aphe- 
lion tangential velocity Vq = 6000 cm s~^. Figure 3 shows the energy 
transfer for a collection of encounters (6^ = 10^ cm) as a function of 

In general, we find that eo > cr, and hence the return orbit of O 
has larger angular momentum than the incoming orbit. While we have 
not rigorously proven this result, it seems intuitively understandable: 
the "lever arm" associated with the Jupiter encounter (i.e., the radius 
of Jupiter's orbit) is ~ 5 times larger than that of the Earth encounter, 
resulting in a larger angular momentum. This finding spoils some of 
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Figure 4- a) Energy transfer per unit mass of O as a function of Earth approach 
distance Be for orbits with Ro = 650 AU and Vb = 6000 cm s~^. b) Eccentricities 
eo (dotted) and cr (sohd) as a function of Earth approach distance bs- 
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the neatness of the scheme. It is possible to reduce the mismatch of 
angular momentum by increasing the Earth encounter distance bE, but 
this change reduces the efficiency of the encounter, as AQe oc 1/bE- 
An example of this behavior is shown in Fig 4. If the return orbit were 
identical to the incoming orbit (modulo its orientation), a mechanism 
could be set up to recycle the object for an indefinite number of passes 
with a very low energy expenditure. 

The angular momentum of O can be restored to its incoming value 
by means of a "course correction" at aphelion. Numerical experimen- 
tation suggests that the most efficient scheme is to adjust bj so as 
to produce a return orbit with the same aphelion as the initial orbit: 
0^(1 + cr) = ao(l + Co). The required velocity change is then simply 
AVr = cr/RO — Vf), applied so as reduce the tangential velocity of O 
to its original value Vq. This velocity correction is similar to the well- 
known Hohmann maneuver used to transfer from one circular orbit to 
another with least velocity change. Since the change AVr is inversely 
proportional to the aphelion distance, it is advantageous (from the 
point of view of least energy expenditure) to arrange for O's orbit 
to have the largest possible aphelion. On the other hand, the orbit 
must not be so large that its period is incompatible with the basic 
encounter timescale of ~ 6000 years which is equivalent to a semimajor 
axis of ~ 330 AU. For typical aphelia 0(600) AU, the velocity change 
is AVr ~ 6000 cm s-^ 

2.2. MULTIPLANET ENCOUNTERS POST-EARTH 

The considerations discussed above suggest that we consider the pos- 
sibility of scheduling multiple planet encounters after passage by the 
Earth. This added complication can help optimize the scheme by re- 
ducing the primary energy expenditure at the return-orbit aphelion. 
An encounter with Saturn, immediately after the Jupiter encounter, 
is a natural candidate. Calculating the post-Saturn orbital parameters 
follows in the manner outlined above. We can then search the parameter 
space of encounter distances with Jupiter bj and Saturn bs to minimize 
the velocity change AVr. 

We find that AVr can in fact be reduced essentially to zero by 
an arrangement where O loses energy at the Saturn encounter; the 
distance bj must be decreased from its best single-encounter value ~ 
1.73 X 10^^ cm to compensate. There is a dramatic decrease of AVr 
to nearly zero (0(10) cm s~^) for bj ~ 6.6 x lO^'^ cm. Figure 5 shows 
the orbital geometry involved. However, to find the minimum AVr 
for any specified value bj demands the specification of bs (or vice- 
versa) to exceedingly high precision. For example, reduction oi AVr 
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Figure 5. Successive encounters (cf. Fig. 2) for an Earth- Jupiter-Saturn arrange- 
ment for an orbit with initial aphelion at 650 AU and aphelion tangential velocity 
Vo = 6000 cm s~^. a) Inner portion of the successive orbits, b) Initial (dotted) and 
return (solid) orbits compared. 



to ~ 10 cm for hj = 6.6 x 10^*^ cm, requires bs to be specified 
to a precision of ~ 10 cm. Less stringent specifications 0(10^) cm are 
sufficient if reduction AVr to a few meters per second is satisfactory. 
Some examples are shown in Fig 6. 

As a result, any realistic scheme will probably not attempt to strictly 
enforce AVr = 0. Nevertheless, it is interesting to find that the "first- 
order" energy expenditure of the scheme can be reduced in principle to 
a negligible amount. Figure 7 shows bj as a function of bg for minimum 
AVr transfers, and also the resulting values of AVr. 



3. Other Issues 

3.1. Accuracy, Course Corrections, Energy Requirements 

In general, various additional effects will interfere with the minimal- 
energy scheme outlined above. These complications include planetary 
orbital eccentricity, non-zero inclination angles, and non-gravitational 
impulses. We will not attempt to develop a sophisticated method of 
guidance for the object, although such methods have been developed 
(to a very high degree of precision) for the space program and plane- 
tary exploration (cf. Battin, 1987). Instead, we will merely make some 
estimates. 

Planetary orbits have eccentricities and inclinations e and i (in ra- 
dians) of a few times 10~^. In order to accommodate these values, the 
velocity changes will be ~ 10^'^V. or about 10*^ cm s~^ when applied 
to velocities of ^ 10 km s~^, which would be typical of the region 
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Figure 6. Velocity AVr required to restore initial orbital energy and angular mo- 
mentum, as a function of difference 6 — hs from Saturn encounter distance for 
three different values of Jupiter encounter distance bj. (Initial aphelion Ro = 
650 AU, aphelion tangential velocity Vo = 6000 cm s~^.) a) bj = 6.8 x 10^°, 
bs = 5.6106641326 x 10^" cm, b) bj = 7.02 x 10^° cm, 6s = 5.8751641993 x lO" cm 
c) 6 J = 7.2 X 10^° cm, bs = 6.1022703373 x 10^° cm 
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Figure 7. a) Jupiter encounter distance hj as a function of Saturn encounter dis- 
tance, required to minimize AVjj. (Initial aphelion Rq = 650 AU, aphelion tangential 
velocity Vb = 6000 cm s"'^.) The upper curve is for encounters in which O gains 
energy at Saturn; the lower for encounters in which O loses energy, b) Resulting 
minimum velocity change AVr for the encounters specified in a). Again, the upper 
curve refers to energy-gaining encounters at Saturn and the lower for energy-losing 
encounters. 



outside of Saturn's orbit. With sufficient planning, one can thus easily 
accommodate the departures of planetary orbits from circles in the 
same plane, as discussed so far. 

Undoubtedly, the need for other course corrections will occur. High 
accuracy is demanded at all critical stages. Closing velocities of 40 km 
s~^ and encounter distances of 10^ cm translate to accuracy of 0(10- 
100) s in time of arrival at the Earth. We can get some feeling for 
the size of velocity changes that are required in a fairly simple way 
by making use of algorithms for the solution of "Lambert's problem" 
(Battin, 1987, Bond and Allman, 1986), which consists of finding the 
velocity vector vq needed to produce a 2-body orbit that takes a mass 
from given position ro to r in a given time interval At. 

We choose a target energy budget for velocity corrections of AF ^ 
10^ cm s"^. We then use the algorithm for Lambert's problem to com- 
pute differential velocity corrections along the incoming orbit (before 
the Earth encounter) that yield changes in arrival times At. We find 
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that a At of 0(10-100) s can be accommodated fairly easily up to 
~ 10^ s before encounter. Alternatively, much larger changes in arrival 
time can be allowed for at greater distances; for example, at ~ 0.5 
AU (~ 4 X 10^ s before encounter), the suggested budget could shift 
the encounter time by ~ 10^ s. Of course, a more realistic mission 
profile would not use up all the allowed energy for one correction, and 
a target velocity change of AV ~ 10"^ cm s~^ per course correction 
might represent a reasonable aim. 

For the case of total AV = 10^ cm s~^, and for a 10^^ gm object, 
about 10^° erg are required to enforce the velocity changes for each 
encounter. About 10^^ erg gm~^ is available from H2O by deuterium- 
tritium fusion, assuming a terrestrial D/H ratio (Pollack and Sagan, 
1993). Thus, ~ 10^^ gm of ice must be proccessed for each encounter (at 
a minimTim); this ice would subtend a volume about 2.2 km on a side. 
For every encounter, this ice volume represents about 10^^ of the mass 
of object O; as a result, the ~ 10^ encounters necessary would consume 
the deuterium of 0{l) large Kuiper Belt object, assuming pure H2O 
ice composition. In addition, about twenty times as much rock mass is 
needed to provide lithium for the production of tritium. An object of 
predominantly chondritic composition may thus be required if in situ 
production is desired. Obviously, much less processing is needed if p-p 
fusion were available; in that case, a single object with an associated 
processing and powerplant could be easily be used for the entire project. 

Non-gravitational forces are also an issue, and they potentially de- 
mand energy expenditures above and beyond those that we have dis- 
cussed thus far. An icy object, although attractive from the standpoint 
of containing fusionable material and favorable initial location in the 
outer solar system, will be more subject to this problem than an stony 
or metallic asteroid. On the other hand, a large main-belt asteroid 
would have to be placed into a suitable orbit starting from the inner 
solar system, where energy requirements are high. 

3.2. Timing 

The successful implementation of this scheme demands a reasonably 
delicate interplay between orbital time scales of a thousand or more 
years (for O) and arrival times scheduled to the minute. The first major 
issue is how often one would expect to have Earth and Jupiter (and 
Saturn) arrive in a particular configuration relative to the argument 
of perihelion of O. The object O spends most of its time "hanging" 
at aphelion. Small adjustments in trajectory can thus be used to time 
the infall to correspond to the moment of proper planetary alignment. 
With this flexibility, we could arrange for O to arrive when Earth and 
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Jupiter are in proper position, an alignment that takes place every 
13 months or so. As a result, two-planct encounters are easily realized; 
three-planet encounters (e.g., including Saturn) are a bit more difficult. 

In order for the most ideal version of our proposed mechanism to 
operate, Earth, Jupiter, and Saturn must all be in the proper phases of 
their orbits when O makes its passage through the inner solar system. 
We can safely assume that minor corrections to the orbit of O can 
delay its arrival in such a way that O's orbital phase is suitable for 
moving the Earth. The three planets, however, must have the proper 
phase relative to each other. The conditions for this phase alignment 
can be written in the form 

{LOq^ — uJ^)t = 27rn , (10) 

and 

{uj(^ — LO']^)t = 2Trk , (11) 

where n and k are integers and where the cuj denote the orbital fre- 
quencies of the planets in obvious notation. 

We first consider the case in which the orbital frequencies, or equiva- 
lently the orbital periods Pj = 27r/a;j, are constant. Expressed in years, 
the planetary periods are approximately P© = 1, Pq^ = 11.86, and Pu 
= 29.28. 

The condition for a perfect alignment can be written in the form 

^J-^Z}}^ (12) 

where we have used the orbital periods rather than the orbital frequen- 
cies. Since the orbital periods are known, the right hand side is a known 
dimensionless number, which has a value of about 18.25. We can write 
this expression in the form 

k 

— = 18.did2dsd4,d5 . . . , (13) 
n 

where the dj denotes digits of the number (which in general will be 
irrational). The general mathematical problem is thus to represent a 
real number (the right hand side above) with a rational approximation. 
For a given specified accuracy (i.e., for a given number of decimal places 
in the above expression), we need a minimum size of the integers k 
and n. The integer k is roughly the number of Earth orbits required to 
attain sufficient alignment, and hence is also the approximately number 
of years between alignments (more precisely, k measures time in units 
of Pe(l - P(b/P%)~^ ~ 1-09 years). 
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One possible choice for the ahgnment integers is thus k = 18did2d3d4^d5 
and n = 10^ , where the last digit represented is the jth one. The time 
interval between alignments is thus about r = P^P'^_n/{PY^ — Pq^) ^ 

20 X 10-? years. 

The above argument shows that a solution exists. A compromise 
must be made, however. In order to increase the accuracy of the align- 
ment, we need longer time intervals between encounters. But we also 
need enough encounters per unit time to move the Earth before the 
Sun compromises the biosphere. 

Although the alignment condition will vary as Earth changes its 
orbital parameters due to the asteroid encounters, the orbital period of 
Earth only changes by a factor of two and hence the accuracy require- 
ments will be of the same order of magnitude for the entire migration 
time interval. The accuracy needed for alignments is determined by the 
accuracy needed for the secondary encounters at Jupiter and Saturn. 
We can assume that the orbit of O will be tuned to interact with 
Earth at just the right impact parameter, but no course adjustments 
can be made before O reaches the outer planets. Jupiter and Saturn 
thus must be in the right place to an accuracy of a few planetary radii 
(say £ ~ 10^*^ cm). This constraint implies that we know the orbital 
phases to a relative accuracy of about £/r ~ 10~^; we must take j = 4 
at the very least, but we would like to use j = 5 or even 6 in an ideal 
case. With j = 4, for example, the time interval between encounters 
(alignment opportunities) is about 200,000 years. During the allowed 
few billion year time period (until the biosphere is compromised), we 
thus only get about 10,000 encounters. But, as discussed above, we 
need nearly a million encounters to successfully move the Earth to a 
viable larger orbit. 

We must thus view the problem the other way around: In order for 
migration to occur within a few billion years, we must have encounters 
every few thousand years. For this frequency of encounters, the largest 
allowed value of the integer k is about 1000 (for example, one obvious 
approximation would be fc=1825 and n=100). With this level of preci- 
sion, we can tune the encounter so that Earth and Jupiter are in the 
right place, but Saturn will generally be in the wrong phase of its orbit 
by an amount corresponding to a fraction ^^0.02 of its orbit. The spatial 
displacement will be 0.02 x27rr ^ 1.2 AU ^ 3000 Rf^ (Saturn planetary 
radii). In this case, we can use Saturn to make course corrections to 
O's orbit , but we cannot obtain perfect post-encounter orbital elements 
(where the object has exactly the same energy and angular momentum 
it started with). 
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Therefore a more realistic goal is to aim at reducing the aphelion 

AV of O to some small value, or at least one that does not dominate 
the "energy budget". There is a range of possible Jupiter encounter 
parameters bj that yield a final AV < 1000 cm s~^; these encounters 
correspond to a range of ~ 0.05 radian of encounters along Saturn's 
orbit, thus easing the timing requirement to a manageable level. 

In addition there are other considerations that mitigate the problem: 

— Uranus and Neptune arc available, giving three times as many 
opportunities as using Saturn alone. 

— Multiple objects can be used for energy transfer, though this will 
probably raise the energy requirements in proportion to the num- 
ber of bodies involved. 

— Encounters need not be scheduled at the first opportunity (as 
shown in Figs. 2 and 5). They can also be timed to occur after 
multiple orbit passes at either intersection point of the orbits. The 
object O can be stored in temporary Chiron-like orbits as well. 



4. Discussion 

In this paper, we have investigated the feasibility of gradually moving 
the Earth to a larger orbital radius in order to escape from the increas- 
ing radiative flux from the Sun. Our initial analysis shows that the gen- 
eral problem of long-term planetary engineering is almost alarmingly 
feasible using technologies that are currently under serious discussion. 
The eventual implementation of such a program, which is moderately 
beyond current technical capabilities, would profoundly extend the time 
over which our biosphere remains viable. 

The main result of this study is a theoretical description of a work- 
able scheme for achieving planetary migration. This scheme is applied 
to the particular case of the Earth. Solar system bodies, such as large as- 
teroids or Kuiper Belt objects, can be used to move Earth over the next 
billion years. These secondary bodies are employed in a gravity-assist 
mechanism to increase the Earth's orbital energy and thereby increase 
its distance from the Sun. The most favorable orbits for the secondary 
bodies have a large semi-major axis, typically hundreds of AU; with this 
relatively high "leverage factor" , the large requisite energy transfer can 
be achieved. 

An important aspect of this scheme is that a single Kuiper Belt 
object or asteroid can be employed for successive encounters. In order 
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to move the Earth at the required rate, approximately one encounter 
every 6000 years (on average) is needed (using objects with mass ~ 
10^2 gm). Due to the acceleration of the Sun's luminosity increase, the 
encounters must be more frequent as the Sun approaches the end of its 
main-sequence life. In order to use the same secondary body for many 
encounters, modest adjustments in its orbit are necessary. However, by 
scheduling the secondary body to encounter additional planets (e.g., 
Jupiter and/or Saturn) in addition to the primary Earth encounter, 
the energy requirements for orbital adjustment at the object's aphelion 
can be substantially reduced. In particular, the energy consumed by 
such course corrections is not likely to dominate the energy budget. 

Any serious proposal for planetary engineering, or any large-scale al- 
teration of the solar system, raises important questions of responsibility 
(sec Pollack and Sagan, 1993). Compared with other astronomical en- 
gineering projects, this scheme has both positive and negative aspects. 
For example, although no massive alteration of planetary environments 
is proposed, this scheme would consume a number of large Kuiper Belt 
objects. 

A great deal of energy must be expended to implement this migra- 
tion scheme. However, the energy needed to move Earth is relatively 
modest compared to that needed for interstellar travel. For example, 
an optimistic minimum energy expenditure is about 10^^ erg, which 
corresponds to the kinetic energy of a ~ 10^^ gm object moving at a 
velocity of 50 km s~^ (this mass is less than 10~^ -^e)- As a means 
of preserving the entire biosphere, this scheme is thus highly efficient 
compared to interstellar migration, even if we have underestimated 
the energy requirements by many orders of magnitude. The energy 
requirements and overall ease of implementation also compare favorably 
with various terraforming projects (Pollack and Sagan, 1993). 

As noted near the beginning of this paper, the required change in 
orbital energy of the Earth is 9 x 10^^ erg. In the basic scheme we 
have outlined, the energy is essentially transferred from Jupiter to the 
Earth. As a result, Jupiter's semi-major axis ai^ decreases by Aa = 
qi^AE/Et^ ~ 5.3 X IQ-^ai^, where E'^ = GMQMn^/2an^ = 1.7 x 10^^ 
erg is Jupiter's orbital energy; this change amounts to ~ 0.03 AU. 
While small, this orbital change could destabilize some asteroidal or 
other orbits by the shift of position of Jupiter's orbital resonances. 
The multi-planet scheme would involve similar-sized orbital changes 
for Jupiter and Saturn (or other planets). 

Potentially more serious questions involve the rotation rate of the 
Earth and the Moon's orbit. We expect that O will raise a tide in the 
Earth during its encounter. The tide could be STibstantial; although O 
would be a relatively small body, the closeness of its passage means 
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that the transient forcing potential would be 0(10) x as strong as that 
of Moon, for a 10^^ gm body passing 10^ cm from the Earth's center. 
Calculating the size and phase of the tide would require detailed work, 
but qualitatively we would expect any tidal bulge to lag in phase behind 
O, as O moves more quickly than the Earth rotates. This in turn implies 
a spin-up of the Earth (similar reasoning accounts for the the spin-down 
of the Earth by the Moon) . Given the very large number of encounters 
planned, a serious increase in the Earth's rotation rate could result. 

However, the above picture, leading to spin-up, takes place only 
for "incoming" encounters, such as depicted in Figs. 1 and 2. The 
symmetry of the encounters equally allows "outgoing" encounters, in 
which O passes by the Earth after its perihelion. Such encounters also 
pass by the Earth's leading limb from inside the Earth's orbit. They 
are thus retrograde with respect to the Earth's rotation, and the same 
considerations as above now lead to spin-down of the Earth rather 
than spin-up. Thus, by careful planning of encounters, we can cancel 
any unbalanced torques exerted on the Earth. 

As for the Moon, reasoning by analogy with cases of stellar binaries 
and third-body encounters suggests that the Moon will tend to become 
unbound by encounters in which O passes inside the Moon's orbit. (As 
well, there is the non-zero probability of collisions between O and the 
Moon, which must be avoided.) Again, detailed quantitative work needs 
to be done, but it seems that the Moon will be lost from Earth orbit 
during this process. On the other hand, a subset of encounters could be 
targeted to "herd" the Moon along with the Earth should that prove 
necessary. It has been suggested (cf. Ward and Brownlee, 2000) that 
the presence of the Moon maintains the Earth's obliquity in a relatively 
narrow band about its present value and is thus necessary to preserve 
the Earth's habitability. Given that the Moon's mass is 1/81 that of the 
Earth, a similarly small increment of the number of encounters should 
be sufficient to keep it in the Earth's environment. 

The fate of Mars in this scenario remains unresolved. By the time 
this migration question becomes urgent. Mars (and perhaps other bod- 
ies in the solar system) may have been altered for habitability, or 
at least become valuable as natural resources. Certainly, the dynam- 
ical consequences of significantly re-arranging the Solar System must 
be evaluated. For example, recent work by Innanen et al. (1998) has 
shown that if the Earth were removed from the Solar System, then 
Venus and Mercury would be destabilized within a relatively short time. 
In addition, the Earth will traverse various secular and mean-motion 
resonances with the other planets as it moves gradually outward. A 
larger flux of encounters might be needed to escort the Earth rapidly 
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through these potential trouble spots. In this case, additional solar 
system objects may require their own migration schemes. 

This technology could also be used, in principle, to move other 
planets and/or moons into more favorable locations within the solar 
system, perhaps even into habitable zones. As an application, the basic 
mechanics of this migration scheme could be employed to clear haz- 
ardous asteroids from near-Earth space. There is also the possibility 
of using Kuiper-belt objects as resources themselves (e.g. of volatile 
materials); gravitational-assist schemes could perhaps deliver materials 
to useful locations with a minimum expenditure of energy. 

An obvious drawback to this proposed scheme is that it is extremely 
risky and hence sufficient safeguards must be implemented. The colli- 
sion of a 100-km diameter object with the Earth at cosmic velocity 
would sterilize the biosphere most effectively, at least to the level of 
bacteria. This danger cannot be overemphasized. 
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